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On  the  Theory  and  Error  Estimation  of  the 
Reduced  Basis  Method  for  Multi-Parameter  Problems^ 

BY 

Werner  C.  Rheinboldt^ 

Abstract.  In  an  earlier  paper  (ZAMM  63,  1983,  21),  J.P.Fink  and  W.C.Rheinboldt  devel¬ 
oped  a  priori  local  error  estimates  for  the  scalar-parameter  case  of  the  reduced  basis  method 
by  considering  the  method  in  a  differential-geometric  setting.  Here  it  is  shown  that  an  anal¬ 
ogous  setting  can  be  used  for  the  analysis  of  the  method  applied  to  problems  with  a  multi¬ 
dimensional  parameter  vector  and  that  this  leads  to  a  corresponding  local  error  theory  also 
in  this  general  case. 


1.  Introduction. 

The  computational  analysis  of  many  problems  in  science  and  engineering  involves  the 
solution  of  very  large  systems  of  nonlinear  equations.  Thus,  not  stirprisingly,  there  is  con¬ 
siderable  interest  in  reducing  the  size  of  these  problems.  One  of  these  reduction  techniques 
hais  become  generziUy  known  as  the  reduced  basis  method.  Some  eztrly  work  on  this  method 
includes  [N77],  [ASB78],  but,  beginning  with  [NP80],  [N81],  [N82],  [NP83a],  much  of  the 
detailed  development  for  a  range  of  problems  in  nonlinear  structinral  mechanics  is  due  to 
A.  Noor  and  his  co-workers.  Since  then  the  literature  on  the  method  has  grown  rapidly 
and  also  the  number  of  application-areas  has  grown  and  includes  now,  for  instance,  initial 
value  problems  for  ordinary  differential  equations  ([PL87]),  differential- algebraic  equations 
([L91]),  certain  classes  of  peu^abolic  differential  equations  ([J90]),  as  well  as  climate  prob¬ 
lems  ([JM871,  (M88]). 

The  method  applies  to  systems  of  parametrized,  nonlinear  equations  and  was  originally 
developed  in  the  setting  of  a  standeird  continuation  procedure  for  nonlinezur  problems  with 
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^Department  of  Mathematics  and  Statistics,  University  of  Pittsburgh,  Pittsburgh,  PA  15260 


1 


one  scalar  parameter.  But,  in  [NP83b],  [NBS84j,  [NPA84],  [NAT87],  A.  Noor  et  al  showed 
that  it  can  also  be  extended  to  problems  involving  several  parameters.  A  first,  genei^l 
analysis  of  the  a  priori  error  behavior  of  the  method  in  the  one-parameter  case  was  given 
in  [FR83].  Some  further  aspects  of  these  results  were  discussed  in  [P85],  and  then  in  [FR85] 
a  general  local  error  estimation  theory  was  developed  which  contains  the  ezurlier  estimates 
as  special  cases.  The  approach  in  [FR83]  was  to  plaice  the  method  in  a  general  differential 
geometric  setting  where  it  can  be  formulated  as  a  projection  onto  a  subspace  spainned  by 
the  first  several  vectors  of  the  moving  frame  of  the  solution  path. 

In  this  paper  we  show  that  this  geometric  approach  can  be  extended  to  the  multi- 
parameter  case.  For  this  we  summarize  in  section  2  the  indicated  differential  geometric 
setting  and  then  show  in  section  3  that  the  general  local  error  theory  of  [FR85]  can  be 
applied  to  the  multi-parameter  case  of  the  reduced  basis  method.  This  is  followed  in  section 
4  with  a  discussion  of  the  basic  approach  to  the  computation  of  the  reduced  basis  vectors 
which  turns  out  to  correspond  to  that  developed  by  A.  Noor  (loc.cit.)  by  different  means. 
At  the  same  time,  the  theory  suggests  some  new  computational  approaches  for  low-order 
approximations. 


2.  The  Basic  Setting. 

-\s  noted  above,  many  stationary  problems  in  techniczd  applications  are  modelled  by  a 
nonlinear  equation 


(2.1) 


F(z,A)  =  0 


f 

'x'  ” 


involving  a  state  variable  z\  a  d-dimensional  parameter  vector  A;  and  a  differential  oper¬ 
ator  F  with  associated  boundary  conditions.  Standard  discretizations  -  for  instance,  by 
finite  elements  or  finite  differences  -  typically  produce  very  large  finite-dimensional  non¬ 
linear  systems  of  equations.  The  aim  of  the  reduced  basis  method  is  to  construct  special 
approximations  that  reduce  the  number  of  degrees  of  freedom  without  compromising  the 
errors. 


□  □ 


The  equation  (2.1)  involves  a  mapping  F  :  Z  x  A  —*  Y  from  the  product  Z  x  A  of  the 
state  space  Z  Jmd  parameter  space  A  into  another  space  Y.  Typically,  the  solution  set 

(2.2)  A/={(r,A)eZxA;  F(z,A)  =  0} 

is  a  differentiable  submanifold  of  Z  x  A.  All  standard  discretizations  leave  the  parameters 
unchanged  and  hence  produce  approximate  equations  of  the  form 

(2.3)  FHizH,\)  =  0 

involving  a  mapping  Fh  '•  Zk  x  A  Yk  where  Zh  C  Z  and  Yh  C  Y  are  stiitably  chosen 
finite  dimensional  subspaces.  The  construction  of  Z^,  Y^,  and  Fh.  usually  talces  little 
or  no  account  of  the  properties  of  the  solution  manifold  M.  Hence,  in  o\ir  differential- 
geometric  terminology,  the  idea  of  the  reduced  basis  method  may  be  viewed  as  the  design 
of  approximations  that  do  take  account  of  information  about  the  manifold  M  and  its  local 
coordinate  systems. 

For  a  classical  illustration  of  this  idea,  suppose  that  M  C  constitutes  a  smooth  curve 

x:  J  C  R^,  X  =  x{s)  E  R^,  sE  J. 

It  is  no  restriction  to  assxune  that  s  denotes  here  the  arclength.  Then  the  local  moving 
frame  at  x  =  i(s)  consists  of  the  tangent  vector  ui(s)  =  x'(s),  the  principal  normal 
U2(s)  =  x"(s)/||x"(s)||2,  and  the  bi-normal  vector  U3(s)  =  ui(s)  x  U2(s).  If  we  express  x 
loczilly  near  x(so)  in  terms  of  this  moving  frame, 

3 

x(s)  =  x(so)  +  ^  Ofc(3,So)ufc(so), 

then  the  coefficients  satisfy  the  well  known  approximation  properties 

aik(s,so)  =  (x(s)  -  x(so))^ujk(so)  =  0((s  -  So)*),  fc  =  1,2,3,  as  s -♦  so- 
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Analogously,  a  suflBciently  smooth  curve  in  il”  can  be  approximated  loczdly  to  order  Jb  + 1 
by  using  as  coordinate  system  the  first  k  moving  frame  vectors,  1  <  it  <  n.  In  other  words, 
independent  of  the  dimension  n  of  the  ambient  space,  edready  a  few  basis  vectors  may 
suffice  to  represent  the  curve  locally  to  a  relatively  high  accuracy.  This  is  the  fundamented 
idea  behind  the  reduced  basis  approximations. 

As  noted  above,  the  given  equation  (2.1)  generally  models  a  problem  in  an  infinite 
dimensional  state  space  Z.  Then,  in  essence,  we  have  the  following  two  possible  choices 
for  introducing  the  desired  approximation: 

(a)  Use  a  ‘reduced-basis’  approximation  of  the  solution  manifold  (2.2)  of  the  original, 
infinite-dimensional  equation  (2.1). 

(b)  Construct  first  a  standard  discretized  equation  (2.3)  of  the  original  problem  and 
then  choose  a  ‘reduced-basis’  approximation  of  the  solution  manifold  of  this  dis¬ 
cretized  equation. 

Evidently,  (a)  has  the  advantage  of  only  one  small  discrete  problem  and  hence  only 
one  discretization  error.  But,  in  general,  this  approximation  is  difficult  to  obtain  since  it 
requires  working  directly  with  the  original  infinite  dimensional  equation.  In  contrast,  (b) 
hzis  the  advantage  that  the  ‘reduced-basis’  approximation  is  much  easier  to  compute  since 
it  involves  only  finite  dimensional  equations.  But  there  is  now  the  disadvantage  that  the 
first  discrete  system  (2.3)  is  large  and  that  the  resulting  two  discretization  errors  are  hard 
to  control  together.  Despite  of  this,  (b)  has  become  the  generally  accepted  approach. 

3.  The  Local  Error  Estimation.  In  line  with  the  approach  (b)  we  assume  from  now 
on  that 

X,  Y  real,  linear  spaces  with  dim  X  =  n,  dim  Y  =  m,  n  =  m  +  d,  d>  1; 

(3.1) 

F  :  S  Y,  C^’-mapping  on  the  open  set  S  C  X,  p>2. 

In  other  words,  F  represents  the  first  discretization  of  the  original  operator.  Generally,  for 
any  linear  space  E  we  denote  an  open  or  closed  neighborhood  of  a  point  e  ^  E  hy  E) 
or  U(e,  E),  respectively.  On  all  linear  subspaces  of  E  the  relative  topology  is  used. 
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As  usual,  a  point  x°  is  a  regular  point  of  the  mapping  F  if  rank  DF{x^)  =  m.  Suppose 
that  as  the  result  of  some  computations  we  have  accepted  a  regular  point  z°  €  5  and 
that  F(x°)  =  y^.  Then  the  submersion  theorem  ensures  that  for  a  sufficiently  small 
neighborhood  U{x°,X)  C  S  the  local  solution  set 

(3.2)  M  =  {  X  €  Liix\X);  F(x)  =  y“  } 
is  a  C^-submanifold  of  X  of  dimension  d. 

In  order  to  introduce  a  local  parametrization  of  this  manifold  at  x°  we  choose  a  splitting 

(3.3)  X  =  T®V,  dim  T  =  d,  dim  r  =  m,  Vn  kerDF(x°)  =  {0}. 

and  consider  the  equation  F(x°  -f-t+u)  =  with  <  €  T  and  v  E  V.  Evidently,  (3.3)  implies 
that  DF{x°)\v  €  Isom  (V,!^).  Hence  the  implicit  frmction  theorem  ensures  the  existence 
of  neighborhoods  U{0,  V)  and  =  U{Q,T)  with  W(x°,-Y)  =  x°  +  +  W(0,  V)  C  S  and 

of  a  C^-map  v  :U^  V  with  t;(0)  =  0  such  that  the  local  coordinate  map 

(3.4)  X,  $(s)  =  x°  +  t  +  u(t),  t  € 

is  a  C^’-diffeomorphism  from  onto  the  (relatively)  open  neighborhood  W(x°,  -Y)  n  M  of 
x°  in  M  (see  also  [R86]). 

A  projection  method  will  be  used  to  construct  the  desired  approximation  of  the  manifold 
M.  For  this  let 

(3.5a)  Pp  €  L{V),  Pj  =  Pp,  rge  Pp  =  Vp,  dim  Vp  =  p  <  m, 

be  a  given  projection  and  note  that,  because  of  dim  V  =  dim  K  =  m  we  can  choose  a 
linear  mapping 

(3.56)  J  6  LiX,  Y),  ker  7  =  T,  rge  7  =  F,  J\v  €  Isom(F,  F). 
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Then  the  subspace  Xp  =  T  ®Vp  C  X  has  dimension  d  +  p,  and  lip  =s  J\vPpJ^y  €  L{Y) 
is  a  projection  of  Y  onto  the  subspace  Yp  =  rge  lip  =  Jjvlp  oi  dimension  p.  Thus  the 
approximate  problem 

(3.6a)  Fp  :  5p  =  5  n  (i°  +  Xp)  ^  Vp,  Fp(x)  =  npF(x),  x  €  5p, 

(3.66)  Fp(x)  =  yj,  x  €  5p,  yj  =  IIpy°. 

is  well  defined. 

Obviously,  x°  €  Sp  not  only  satisfies  the  original  equation  F(x)  =  y°  but  also  the 
approximate  equation  (3.6b).  We  are  interested  in  the  solutions  of  (3.6b)  locally  near  x° 
and  in  am  estimate  of  the  approximation  error.  For  this  we  extend  the  mapping  (3.6a)  to 
eJl  of  5  in  line  with  the  following  result: 

Lemma  1.  Consider  the  extended  mapping 

(3.7)  Fp  :  5  y,  Fp(x)  =  (/y  -  np)J|vr(x  -  x“)  +  np(F(x)  -  y°), 

where  F  €  L(X)  denotes  the  natural  projection  of  X  =  T  0  F  onto  the  second  factor. 
Then  Fp(x)  =  0  for  some  x  €  5  if  and  only  if  x  €  Sp  and  Fp{x)  =  y°. 

The  proof  is  straightforward.  In  fact,  if  Fp(x)  =  0forx  =  x°+t  +  y  G  5  then 
J|vy  =  llp{J\yv  —  F(x)  +  y°)  €  Yp  whence  v  ^  Vp  eind  therefore  x  €  Sp  as  well  as 
( Jy  —  rip) J|yu  =  0  and  Fp(x)  =  y®.  Conversely,  if  Fp(x)  =  y®  for  i  =  x°  + 1  + 1?  €  Sp  then 
V  G  Vp  and  thus  Fp{x)  =  0. 

Evidently,  we  have  Fp{x°)  =  0  and  it  turns  out  that  x°  is  a  regulzir  point  of  Fp  if  the 
condition 

(3.8)  Fp  n  ker  np£>F(x°)  =  {0}, 
holds.  In  fact,  by  (3.7)  we  have 

DFpix°)  =  {ly  -  np)j,yr  +  npF>F(x°) 
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and  (3.8)  shows  that  DFp{x°)\v^  €  Isom  {Vp,  Yp)  whence  IIpI? F(z°)  maps  V  onto  Yp  while 
(7j,  —  Up)  Jiv^r  maps  V  onto  the  complement  rge  (/p  —  Up).  Therefore,  all  points  z  in  some 
sufficiently  small  neighborhood  i/(z°,-Y)  C  S  of  z°  are  regular  £ind  the  local  solution  set 

(3.9)  Mp  =  {  z  €  U{x\X);  Fp(z)  =  0  }, 

is  again  a  d-dimensional  submanifold  of  X  of  class  C’. 

It  turns  out  that,  because  of  (3.8),  we  can  define  the  approximation  error  by  comparing 
points  on  M  and  Mp  near  z“  that  have  the  same  coordinate  t.  This  leads  to  the  following 
theorem  proved  in  [FR85]  (see  also  [R86]): 

Theorem  1.  Under  the  assumptions  (3.1)  let  z°  €  5,  with  F(z“)  =  y°,  be  a  regular  point 
where  on  the  local  manifold  M  of  (3.2)  the  splitting  (3.3)  dehnes  the  local  coordinate  map 
(3.4).  With  the  projection  (3.5a)  and  the  linear  map  (3.5b)  introduce  the  approximate 
problem  (3.7a/b)  and  suppose  that  (3.8)  holds.  Then  there  exists  a  ball  Bq  C  centered 
at  the  origin  ofT  and  a  neighborhood  14(0,  V)  such  that  x°  +Bo  +14(0,  V)  C  S  and  that  for 
any  t  ^  Bo  the  equation  Fp(x°  +  t  +  w)  =  0  has  a  unique  solution  w  =  Vp(t)  in  14(0,  F).  The 
mapping  Vp  :  Bo  — ►  W(0,  V)  is  of  class  and  satisfies  Up(Bo)  €  Vp  ni4(0,V).  Moreover, 
the  error  estimate 

(3.11)  ||i;(t)  -  up(<)l|  <  C\\(Iv  -  Fp)y(OII,  V<  e  Bo, 

holds  with  a  constant  C  independent  of  t. 

We  sketch  briefly  the  proof  zuid  refer  for  further  details  to  [FR85]  or  [R86].  It  is  based 
on  the  following  implicit  function  theorem  used,  in  a  more  general  form,  in  [BRR81a,b,c]: 

Theorem  2.  Let  V,  T,  1*  be  Snite  dimensional,  real,  linear  spaces  and  G  :  S  Y  a  C^- 
p  >2,  on  an  open  set  S  of  X  =  T  xV  with  (0,0)  €  S  and  G(0,0)  =  0.  Suppose  that 
on  a  compact  neighborhood  Uc  =  U(0,T)  there  exists  a  Lipschitz  -  continuous  function 
V  :  Uc  —*  V  such  that  v(0)  =  0.  t  +  v(t)  G  S  for  t  €  Uc,  and 
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(i)  D„Git,v{t))€lsomiV,Y),^t£Uc, 

(ii)  ||I>„G(<,t;(0)-‘|l  <  71,  ||A(?(#,«(0)ll  <  72,  V<  6  Uc, 

(Hi)  I1G(<,  u(<))ll  =  0, 

where  71,  72  are  positive  constants  independent  of  t.  Then  there  exists  a  ball  Bo  C  Uc 
centered  at  the  origin  of  T  and  a  neighborhood  U(0,  V)  such  that  Bq  ®  W(0,  V)  C  5  and 
that  for  any  t  €  Bq  the  equation  G(t,w)  =  0  has  a  unique  solution  w  =  w(t)  €  U(0,V). 
Moreover,  the  mapping  w  •.  Bo  —*V  is  of  class  C*  and 

(3.12)  ||u;(t)-i;(0)||  <C||t||,  Vt  €  Bo 

with  a  constant  C  independent  of  t. 

The  main  part  of  the  proof  of  Theorem  1  consists  in  applying  Theorem  2  to  the  mapping 

G:Sg=  t€T  xV,x°  +  t  +  v€  5}  -*Y,  G{t,v)  =  Fp(x°  +t  +  v). 

By  Lemma  1  we  have  G{t,v)  =  0  for  (t,v)  €  So  exactly  if  i  =  +  t  +  u  6  5p  and 

Fp(i)  —  1/p.  Moreover,  (3.8)  implies  that  D„G(0,0)  €  Isom  {V,Y).  Hence 

D„G{t,v)^  Isom  (YY),  ||£>„G(t,  i;)-'||  <ai,  ||I>tG(t,  i;)||  <  oj, 

for  adl  (t,  u)  in  some  compact  neighborhood  F{(0, 0),  T  x  V)  C  Sc-  We  restrict  this 
neighborhood  such  that  for  all  (t,  ui ),  (t,  V2)  €  I7((0, 0),r  x  V') 

||G(<,U2)  -  G(t,vi)  -  D„Git, vi){v2  -  i;i)||  <  t^||u2  -  ui|| 

20, 


whence  it  follows  readily  that 

(3.13)  ||t;2  -  II  <  Y^lGit,  V2)  -  G(t,  )|| ,  V(t,  y, ),  (t,  )  €  t/((0, 0),  T  x  V), 
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Since  the  mapping  v  :  —*  V  in.  (3.4)  is  of  class  C*  and  u(0)  =  0  there  exists  a 

compact  neighborhood  Uc  =  U{0,  T)  C  such  that  (f,  v{t))  €  t^((0,0),  T  x  F)  for  t  €  Ue- 
This  imphes  that  the  conditions  (i)  and  (ii)  of  Theorem  2  hold  on  Uc-  Moreover,  (iii)  is  a 
direct  consequence  of 

=  ll(/y  -  n,)J|vr(<K()  -  x”)!!  <  const  ||($(()  -  i'’)||. 

Thus  Theorem  2  applies  and  there  exists  a  ball  Bq  C.  centered  at  the  origin  of  T  and  a 
neighborhood  U(0,  V)  such  that  +Bo  +U(0,  V)  C  S  and  that  for  any  t  €  Bq  the  equation 
Fp(x°  +  t  +  lu)  =  0  has  a  unique  solution  w  =  Vp{t)  in  U(0,  V).  Moreover,  the  mapping 
Vp  :  Bq  —*  U{Q,  V)  is  of  class  and  it  follows  from  Lemma  1  that  Up(Bo)  €  Vp  nl/(0,  V). 
Finally,  (3.12)  ensures  that  ||up(t)||  <  C'||<||  for  t  €  Bq.  Hence,  by  restricting,  if  needed,  the 
ball  Bq  we  have  €  C/'((0,0),r  x  F)  and  the  estimate  (3.11)  follows  directly 

from  (3.13). 

4.  The  Choice  of  the  Reduced  Basis  Vectors. 

The  local  error  estimate  (3.11)  suggests  that  we  should  construct  the  approximating 
space  Fp  such  that  the  quantity 

(4.1)  ||(/v  -  np)t;(#)|l 

is  as  small  as  possible.  For  this  we  consider  at  the  given  regular  point  S  the  ‘tzmgential’ 
splitting 

(4.2)  .Y  =  T©F,  T  =  ker  I>F(i°),  F  =  rge  I>F(i”)^. 
for  which  (3.2)  certainly  holds.  Here  we  may  use 

(4.3)  7  =  DF(i°)€  KXF). 
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as  the  mapping  (3.5b).  The  term  ’tangential’  splitting  derives  from  the  fact  that  the 
tangent  space  of  the  local  solution  manifold  (3.2)  of  F  at  x°  may  be  identified  with 
ker  DF{x°). 

For  the  computation  we  need  basis  mapping  s 

(4.4)  A'€l(i2‘',A'),  rgeA'  =  T, 
and 

(4.5)  A  €  X(ii’’, *Y),  rank  A  =  p,  rge  A  =  Vp  CV. 

of  the  local  parameter  space  T  and  the  approximating  space  Vp,  respectively.  By  the  choice 
(4.3)  of  J  we  have  JA  €  L(R’’,Y)  and  rank  JA  =  p  whence  there  exists  a  bi-orthogonal 
mappping 

(4.6)  B€L(RP,YI  B^JA  =  Ip, 
eind,  the  corresponding  projections  Pp  eind  Up  can  be  written  as 

(4.7)  Pp  =  AB'^J,  Tip  =  JAB'^. 

Therefore  the  resulting  approximate  equation  is  equivalent  with  the  reduced  problem 

(4.8)  G{s,  w)  =  B^F(x®  +  Ks  +  Aw)  =  b°  =  B^y°,  G  :  R^  x  R’^  A". 

The  condition  (3.8)  is  here  automatically  satisfied.  In  fact,  if  TlpDF{x°)u  =  0  for  some 
u  ^  Vp  then,  since  JAw  =  0  only  if  iw  =  0,  we  have  B^DF{x°)u  =  0  and  therefore 
ker  DF(x°)nVpC  ker  ^^(x")  H  V  =  {0}. 

For  the  representation  of  the  local  coordinate  map  (3.4)  it  is  useful  to  introduce  the 
mapping 

(4.9)  77  :7/(0,A'') ->  r,  77(5)  =  £>F(x'’)+u(As),  s  €  W(0,  A**), 
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where  as  usual  DF{x°)'^  =  {DF(x°)DF{x^)^)~^ DF{x°)  and  hence  DF{x°)DF{x'^)'^  is 
the  orthogonal  projection  onto  V.  Then  we  have  v(Ks)  =  DF{x^)^t}{s)  for  s  €  U{0,R'^). 
The  tangency  property  of  the  splitting  obviously  implies  that  Dt}{0)  =  0.  Thus  for  any  k 
less  thsin  the  differentiablity  class  p  of  F  v/e  obtain 

*  1 

(4.10)  viKs)  =  Y,  -DFixYD^niO)s^^^  +  ^(^),  ^  €  Uo  =U(0,R'^) 

j=2^' 

where  the  remainder  R(s)  satisfies 

(4.11)  lIRWII  =0{||^||‘'4‘),  as||»||-.0 

and  we  used  the  notation  =  (s,s,...  ,s),  (j-times). 

Therefore,  if  we  choose  the  approximating  space  Vp  such  that  for  some  k  >2 

(4.12)  DFixYD^v(0)s^^^  €Vp,  >  =  2, . . .  ,  ib,  Vs  € 
then  it  follows  from  (4.10)  and  (4.11)  that 

(4.13)  ||(/v,-np)i,(A's)|l=l|R(3)|l=0(||3||‘4‘),  asl|3l|^0, 

and  we  obtain  from  Theorem  1  the  desired  asymptotic  estimate 

(4.14)  ||i-(A-»)  -  „p(A-3)||  =  0(||s||*4').  as  ||s||  ^  0. 


For  fixed  s  €.  Uq  the  mapping 

(4.15)  r  €  R^,  \t\<6  ^  7,(r)  =  x®  +  tKs  +  DF{x°fTf(Ts) 

defines  a  path  on  M  through  x°  which  has  at  x°  the  derivatives 

(4.16)  D7,(0)  =  t  =  Ks,  D^jAO)  =  DF{x^)D^r,{Q),  j  =  2,...,  p. 


11 


The  natural  inner-product  of  X  induces  on  Mo  a  Riemannian  structure  zmd  the  derivatives 

(4.16)  reflect  geometric  properties  of  this  Riemannizui  manifold  near  z°.  For  example,  the 
second  orde,.  term  satisfies 

(4.17)  =  -DF(i°)^[DF(i°)DF(i“)^]-»I>2F(i°)t(2)  =  v(i,t),  t  =  Ks^T 

where  V  denotes  the  second  fimdamental  tensor  of  Mq  at  (see  [RR89])  and  hence  this 
term  reflects  cmrvature  properties  of  Mq.  Similarly 

DF(x°fD^ri(0)s^^^  =  -f-  3D^F(x°)(t,  V(t,t)  -  DF(x°)t, 

represents  a  first  ’’torsion”  tensor. 

Since  F('Yj(t))  =  0  for  sufficiently  smaill  r  6  it  follows  by  repeated  differentiation 
that 

(4.18)  DF(x°)£f^7,(0)s^^^  =  ff;(x'’)<(»,  f  =  A's  €  T,  ;  =  2, . . .  ,  it, 

with  multi-linear  operators  Hj(x^)  that  are  readily  calculated  recursively.  For  example, 
we  have 


H3{x^)t^^^  =  -D^F(x‘’)t<3>  -|-3I>2F(i°)(t,i?S,(0)s(2)). 


For  the  computation  of  the  derivatives  occurring  in  note  that  be¬ 

cause  of  K'^7,{t)  —  K^(x°  +  tKs)  for  r  €  iZ**  we  have  A' ^F-'7,(0)  =  0  for  j  =  2, . . .  ,  it. 
Therefore,  these  derivatives  tire  the  unique  solutions  of  the  systems 


(4.19) 


rn 


t  =  Ks  €  T,  j  =  2,...  ,Jk 
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since,  by  construction  of  A',  the  matrix 

is  nonsinguiar.  The  systems  (4.19)  were  derived  differently  by  A.Noor  et  aJ  (see  e.g. 
[NR83])  who  also  showed  them  to  be  an  effective  tool  for  the  computation  of  the  needed 
path  derivatives  for  many  practical  problems.  For  a  related  discussion  about  some  of  the 
numerical  aspects  we  refer  also  to  [M87]. 

For  the  construction  of  the  approximating  subspace  Vp  C  X  that  satisfies  the  condition 
(4.12),  let  si ,  S2?  •  •  •  be  some  basis  of  and  compute  for  given  k  the 


(4.20) 


P  = 


-d-1 


vectors 


(4.21)  +  =  2, . . .  ,  fc. 

These  vectors  span  the  desired  subspace  Vp  C  X.  In  practice,  we  usually  finds  that  its 
dimension  p  equals  p  in  which  case  the  vectors  (4.21)  can  be  used  directly  as  the  columns 
of  the  matrix  A.  Otherwise,  there  axe  various  means  of  extracting  a  linearly  independent 
subset  to  construct  A. 

Besides  the  p  vectors  (4.21)  we  need  d  vectors  for  AT.  In  view  of  (4.20)  the  total  number 
d  +  p  of  vectors  that  have  to  be  computed,  increases  very  rapidly  as  the  following  table 
shows:  _ _ _ _ _ 


k 

d=  1 

d  =  2 

d  =  4 

d  =  6 

d=  10 

1 

1 

2 

4 

6 

10 

2 

2 

5 

14 

27 

65 

3 

3 

9 

34 

83 

285 

4 

4 

14 

69 

209 

1000 

For  przictical  applications  this  limits  the  reduced  basis  approach  to  problems  with  a 
smzdl  number  d  of  parameters  and  to  a  choice  of  small  orders  of  approximation  k.  But  the 
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use  of  small  fc,  say  Ar  =  2  or  fc  =  3,  obviously  requires  more  frequent  updates  of  the  bzisis 
vectors  which  in  turn  calls  for  the  development  of  faster  algorithms  for  the  computation 
of  these  vectors. 

We  conclude  this  discussion  by  showing  that  for  k  =  2  the  relation  (4.6)  with  the  second 
fundamental  tensor  of  M  leads  to  a  direct  algorithm  for  computing  K,  A,  and  B  based 
on  the  results  of  [RR90].  For  this  suppose  that  natural  bases  are  used  on  X  =  and 
y*  =  iZ"*.  With  the  QR- factorization 


where  rge  Q2  =  ker  DF(x°)^,  we  can  set  K  =  Q2-  For  any  given  s  £  with  sufficiently 
small  ||<||  <  e  let  X  =  +  Ks.  Then  the  chord  Newton  algorithm 

while  ”no  convergence” 

solve  RT z  =  F{x)  for  2; 

set  X  :=  X  —  Q(r,0)^; 

converges  to  $(s)  =  +  Ks  +  DF(x’^)^Tf{s)  €  M  and  hence  implements  the  tzmgential 

coordinate  system. 


The  algorithm  of  [RR90]  for  computing  the  second  fimdamental  tensor  aissiimes  that 
the  natural  inner  product  is  used  to  induce  a  Riemanniain  structure  on  the  mamifold  M. 
Briefly,  for  d  =  2  and  fc  =  2  let  ui  and  U2  denote  the  two  columns  of  K  =  Q2  and  set 
U3  =  ui  4-  U2.  It  was  shown  in  [RR90]  that  for  u  =  Ui  with  fixed  index  i  =  1, 2, 3,  and  any 
sufficiently  small  /i  >  0  the  algorithm 


yi  =  ^{h) 


y2  := 

V  :=  (yi  -  x°)  +  {y2  -  x°) 

V  :=  V  —  ujvui 

V  :=  t;/||v|| 


chord  Newton  process 
chord  Newton  process 
approximation  of  the  normal  vector 
orthogonalization  w.  r.  t.  Ui 
normalization  to  length  1 


produces  an  approximation  of  u,-  =  V{ui,ui).  With  W  =  (ui,  ^2,^3)  we  compute  now  the 
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singular  value  decomposition  of  the  3x3  matrix  JW: 


Bl 

BJ. 


JW{C^,Cx)  = 


diag  ,<Tp)  o' 
0  0 


Here  we  have  Bq  €  and  Co  €  and  p  <  3  denotes  the  maximal  index  such  that 

(7i  >  . . .  >  <Tp  >  0.  Then  it  follows  that  p  =  rank  A  and  hence  that 


A  =  WCq^  5  =  Bo  diag  ,l/<Tp). 

are  the  desired  basis  matrices. 

The  algorithm  can  be  efficiently  integrated  into  the  method  for  computing  simplicial 
approximations  of  M  developed  in  [R88].  This  also  opens  up  the  possibility  of  extending 
the  approach  to  the  case  d  =  2,  =  3  by  using  finite  difference  approximations  for 

computing  the  required  third  order  basis  vectors.  The  resulting  overall  process  will  be 
explored  elsewhere. 
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